Dirac Fermion Confinement in Graphene 
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We study the problem of Dirac fermion confinement in grapliene in the presence of a perpendicular magnetic 
field B. We show, analytically and numerically, that confinement leads to anomalies in the electronic spectrum 
and to a magnetic field dependent crossover from \/B, characteristic of Dirac-Landau level behavior, to linear in 
B behavior, characteristic of confinement. This crossover occurs when the radius of the Landau level becomes 
of the order of the width of the system. As a result, we show that the Shubnikov-de Haas oscillations also change 
as a function of field, and lead to a singular Landau plot. We show that our theory is in excellent agreement with 
the experimental data. 
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The production of two-dimensional (2D) graphene 011^, 
and the confirmation, via an anomalous integer quantum Hall 
effect 1 3, 4], of the presence of Dirac particles in its electronic 
spectrum, has attracted a great deal of interest. Because of 
the vanishing of the density of states at the Dirac point, these 
semi-metallic systems present properties that deviate consid- 
erably from Landau's Fermi liquid theory j^^. In fact, these 
systems show properties that are similar to models in parti- 
cle physics and, in particular, to relativistic quantum electro- 
dynamics (QED) but with an effective "speed of light" (the 
Fermi-Dirac velocity, vp) that is substantially smaller than the 
actual speed of light, c (wf ~ lO^^c). In the most general 
case, the electron dispersion in graphene can be written in the 
form of Einstein's equation: E±{k) — ±y/m'^Vp + Vpk'^, 
where k is the electron momentum (from now on we use units 
such that c — h = 1 = fc^) and m is the relativistic mass. 
In solids this mass represents a gap, A = mvp, in the elec- 
tronic spectrum. This gap can be generated, for instance, by 
the spin-orbit coupling tJ . 

Furthermore, due to experimental constraints, graphene 
samples are usually mesoscopic in size |8, 9] leading to a sit- 
uation where Dirac fermions are confined by either zig-zag 
or armchair edges to a finite region in space 1 10|. Confine- 
ment is also particularly important for the production of elec- 
tron wave-guides that are the main elements for the produc- 
tion of electronic devices such as all-carbon transistors. Dirac 
fermion confinement was a particularly enigmatic problem in 
the early days of quantum mechanics since the formation of 
wave packets in a region of the size of the Compton wave- 
length, w 1/m, requires the use of negative energy solutions, 
or anti-particles, leading to a ground state with time dependent 
cuiTents, the phenomenon called zitterbewegung 111 ill . An- 
other manifestation of this confinement effect is Klein's para- 
dox where a flux of particles incident on a square potential 
baiTier produces a reflected current that is larger than the inci- 
dent one. 

In this paper we show that graphene's zitterbewegung can 
be studied directly with the application of a transverse mag- 
netic field. We show that the confinement, generated by the 



finite size of the sample, shows up in a rather non-trivial way 
in the electronic spectrum. In particular, we show that the so- 
called Landau plots (the dependence of electronic spectrum 
in the magnetic field |12)) is rather non-trivial when the cy- 
clotron length becomes of the order of the size of the sample. 
We address this problem analytically by studying the Dirac 
equation in a magnetic field and also by solving numerically 
the tight-binding model for graphene in a finite geometry. 

Graphene is a honeycomb lattice of carbon atoms (with two 
sublattices, A and B) with one electron per 7r-orbital (half- 
filled band) and can be described by a tight-binding Hamilto- 
nian of the form: 



t.b. 



(1) 



where a\ ^ (ai^a) creates (annihilates) an electron on site R^, 

with spin a {a =t,i), on sub-lattice A, bl ^ (bi^a) creates 
(annihilates) an electron on site R^, with spin a, on sub-lattice 
B, tij = tex-p{i9ij} is the nearest neighbor hopping energy 
(t w 2.7 eV) in the presence of a magnetic field H — B z 
(9ij = 2tt A ■ dl/^o, with A = Bxy and $o = 2n/e is 
the quantum of magnetic flux). In the absence of next-nearest 
neighbor hopping, t' (w O.lt), the Hamiltonian is particle- 
hole symmetric |6] (the Zeeman energy is disregarded). 

In a finite system, one has to add the confining poten- 
tial: He = J2i 1 where rii is the local electronic den- 
sity. Vi vanishes in the bulk but becomes large at the edge 
of the sample. We have studied different types of potentials 
(hard wall, exponential, and parabolic (3) but in this paper 
we will focus on a potential that decays exponentially away 
from the edges into the bulk with a penetration depth, A. In 
Fig. n we show the electronic spectrum for a graphene rib- 
bon of width L = 600a (a — 1A2 A, is the Carbon-Carbon 
distance), in the presence of a confining potential, V{x) = 
Vo [e-(^-'f'/2)/^ + e-('^/2-^)/^], with strength Vq = O.lt, 
and a penetration depth A = 150a (we choose this large value 
of A just to illustrate the effect of the confining potential in 
detail, in real samples we expect A w a, which is the case dis- 
cussed in the text), as a function of the momentum k along the 
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ribbon. One can clearly see that in the presence of the con- 
fining potential the particle-hole symmetry is broken and, for 
Vi > 0, the hole part of the spectrum is strongly distorted. In 
particular, for k close to the Dirac point, we see that the hole 
dispersion is given by: En,a=-i{k) ~ —Jnk'^ — Pnk'^ where 
n is a positive integer, and 7„ < (7„ > 0) for n < N* 
(n > N*). Hence, at n = N* the hole effective mass diverges 
(jN* = 0) and, by tuning the chemical potential, via a back 
gate, to the hole region of the spectrum (/i < 0) one should be 
able to observe an anomaly in the Shubnikov-de Haas (SdH) 
magneto-transport oscillations. 
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FIG. 1 : (color on line) Energy spectrum (in units of t) for a graphene 
ribbon 600a wide, as a function of the momentum k along the rib- 
bon (in units of l/(\/3a)), with a magnetic flux of 5 x lO^^^o per 
hexagon, in the presence of confining potential (see text). 

At low energies and long wavelengths, the energy spectrum 
of Hamiltonian ([l) reduces to two Dirac cones centered at the 
K and K' points in the Brillouin zone. Around each Dirac 
point the Hamiltonian Q can be written as: 



Hq = a- ■ (vfP + eA) . 



(2) 



where vp = 3ta/2, or are Pauli matrices acting on the 
states, ■(/'(r) = ('!/'A(r), '0_B(r)), of the two sub-lattices, and 
p = {px,Py) = — iV is the 2D momentum operator In the 
absence of confinement (V{x) = 0) we can diagonalize (|2j 
and one finds Landau levels given by: 



En, a = (J\pivp\pn jlu — crvFV^^eB n , 



(3) 



where £b = l/veB is the cyclotron length, n is a positive 
integer, and a = 1 (—1) labels the electron (hole) levels. 

As discussed in the context of neutrino billiards 1 14], the 
problem in the continuum suffers from the difficulty that in 
trying to confine massless Dirac particles in a region of size 
L by including a large potential V at the edge, leads to a sit- 
uation where particles still exist even at energies higher than 



V. This problem, of course, does not arise in the tight-binding 
description. In order to avoid this problem in the continuum 
description, we introduce a position dependent mass term: 
He = VpM{x)<7z , where, 



Mix) - <^ 0, 



M, x<~L/2, 
<x< L/2, 
M, x>L/2, 



(4) 



where L is the width of the graphene stripe. We are interested 
in the hard wall case (M oo) although other potentials can 
be studied in analogous way 1 13]. Notice that in the absence of 
an applied magnetic field, a mass term does not break particle- 
hole symmetry, as in the case of a potential F(a:))|[3l- Never- 
theless, since both V{x) and M (x) are strongly concentrated 
at the edge (in a distance of the order of the lattice spacing), 
they do not modify the states in the bulk. It is also worth men- 
tioning that although (|3 is not time reversal symmetric in the 
absence of a magnetic field [14], time reversal is recovered by 
the inclusion of the second Dirac cone at the opposite side of 
the Brillouin zone. 

The Dirac equation. Hip — Etp, where H = Hq + He, can 
be recast in terms of a wavefunction ansatz, <j) = (H + E)ip, 
as: H^<j) — E'^cj). It is easy to show that this wavefunction has 
the form: y) cx e'^''^^ ipcr{x) , where k is the momentum 
along the y direction, cr = ±1 are the eigenstates of (Tz, and 
(pcr [x) obeys the following equation: 



where. 



v{0 



{iBE^jvpf - 



(5) 



(6) 
(7) 
(8) 



Equation (|5} is a dimensionless Schrodinger equation for a 
non-relativistic particle (of mass 1/2) in a parabolic potential 
(of frequency Wo = 2) superimposed to a potential well, y(^), 
whose position shifts with the momentum k (see Fig.|2j- We 
see, from (|8j, that the Dirac fermion spectrum in the presence 
of the magnetic field and confining potential can be written as: 



(9) 



At low energies (small e) the parabolic potential dominates 
and the wavefunctions look like the ones in the infinite sys- 
tem {L oo) in the presence of magnetic field |6]. In this 
limit, the spectrum of (|5} is given by the ID harmonic oscil- 
lator: e ^ u!Q{n + 1/2) = 2n + 1. This result, together with 
(|9j, gives rise to eq. Q, with the energy proportional to ^/B. 
On the other hand, at large energies, the confining potential 
becomes more important and the energy spectrum changes to: 
e « [7rn/(i/^s)]^, leading to a Dirac fermion spectrum of 
the form: 



E„ 



avF[T^n/ L + aLeB /{2Tm)] 



(10) 
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FIG. 2: (color on line) Illustration of the potentials in the problem: 
= (top) and fc / (bottom). 

which leads to a spectrum proportional to B. These simple 
arguments show that the SdH magneto-resistance oscillations 
changes behavior as a function of magnetic field. On the one 
hand, for a given chemical potential /i, eq. (|3ji predicts that the 
maxima of the SdH should happen at fields: 

l/BF{N,L)^{2vle/t?)N , (11) 

where N is the Landau level index. On the other hand, eq. (I10> 
shows that the maxima occur at fields: 

1/Bf{N, L) = [L'e/{2^''N)]/{N,{L) - N) , (12) 

where 

NciL) - ^iL/invp) . (13) 

Hence, l/Bp{N) diverges at a critical Landau level index 
Nc{L) which increases linearly with the width L of the 
graphene stripe. The deviation from (II 1> to (I12> is a clear 
sign of the Dirac fermion confinement. 

Notice that the crossover from Q to (I10> (or from (II 1> and 
(I12» occurs when the Landau orbit fits into the confining po- 
tential. Since each orbit must enclose exactly an integer num- 
ber of flux quantum, $o> the crossover occurs at a magnetic 
field B* such that, B*L'^ w A^^o, that is, 

B*{N,L) N^o/L^ . (14) 

Let us now consider the numerical solution of the differen- 
tial equation (|5j written in terms of the above introduced di- 
mensionless variables in the case M oo. Because our treat- 
ment of the Dirac equation leads to a second order differential 
equation, the appropriate boundary conditions for a sharp con- 
fining edge with a mass term is 1 15 1: (pa{ — L/{2£B) + k£B) = 
<fcr{L/ {21b) + fc^s) ~ 0. In Fig.|3lwe show the energy spec- 
trum at i? = 1 T for two different system sizes as a function of 
kP' jL. One can clearly see that the degeneracy of the Landau 



FIG. 3: (color on line) Energy spectrum, Ecr{k), as function of 
k£%/L for two system sizes, L = 50 nm and L — 500 nm, at 
B = 1 T, for 10 eigenstates. Each Landau level is sub-divided into 
two sub-levels, distinguished be the quantum number a. 
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FIG. 4: (color on line) Wave functions ipa for L = 500 nm. The 
values of the magnetic field are, from left to right and top to bottom, 
5 = 5, 2.5, 1.25, 0.6, 0.3, 0.1, 0.05, 0.025, 0.01 (T). 

levels is lifted for small enough system sizes or large enough 
k (l^ UtIi . For small k the energy states are dispersionless 
and degenerate. 

In Fig.|4]we show the first two state eigenvalues of the effec- 
tive Schrodinger equation (|5}, for fc = and different values 
of the field (or, equivalently, different system sizes). One can 
clearly observe the change in the wavefunction from cosine 
(sine) to gaussian (first order Hermite polynomial times gaus- 
sian) behavior as Ib decreases. Clearly the system evolves 
from a state where the boundaries introduced by the confin- 
ing potential are iiTelevant (the wave functions and the energy 
levels are essentially those of the ID harmonic oscillator at 
B = 5 T, and /s ~ 12 nm), passing to a state where the gaus- 
sian decay of the wave function in the classically forbidden 
regions is important allowing the electrons to feel the pres- 
ence of the confinement potential (the wave functions and the 
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energy levels cannot be described either by the ID harmonic 
oscillator or by the particle in a box for i? = 0. 1 T at ~ 80 
nm). Finally, when the Landau orbit is of the order of size of 
the confinement potential, the eigenstates are essentially those 
of the particle in a box {B = 0.01 T, Is — 250 nm). 
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■ Experimental data 
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FIG. 6: (color on line) On the left hand side, we show the energy 
spectrum (in units of t) as a function of (in T~^). The hori- 
zontal line shows position of the Fermi energy (~ 0.069t). On the 
right hand side, we show the theoretical Landau plot (open circles) 
in comparison with the experiments of ref. L18.1 (close squares). 



FIG. 5: (color on line) Landau plots for two system sizes, L= 500, 
50 nm. On the right hand side, we depict the energy levels (in K) as 
a function of magnetic field (in T) and the horizontal line marks the 
position of the Fermi energy. On the left hand size we plot 1 / Bf{T) 
as a function of the Landau index N . 

In Fig.|5] we show the energy spectrum as a function of the 
magnetic field for different system sizes together with their 
respective Landau plots. Note that at small fields (when L 
is large enough) the energy spectrum follows the ^/B depen- 
dence of (|3} while at larger fields it becomes linear in B as 
predicted by The crossover from these two asymptotic 
behaviors is indeed given by eq. \\A\ . as one can see from 
the size dependence. More striking, however, is that fact that 
\/Bf indeed diverges at sufficiently high Landau level index 
and that the size dependence is given by ( I13> . 

In Fig.|6lwe compare our tight binding results with with the 
experimental data of ref. 118]. We choose a ribbon of size 295 
nm (equivalent to A^c — 1197 unit cells), and Fermi energy 
0.069t (equivalent to 0.22 eV). Notice the excellent agreement 
between theory and experiment for < 16 (2T < B < 
lOT). For > 16 there is shift of N by one (either plus 
or minus one) relative to the experiment. This discrepancy, 
we believe, can be assigned to the experimental difficulty in 
assigning the Landau indices at small magnetic fields 1 18|. 

In summary, we have studied the problem of Dirac con- 
finement in graphene, that is, graphene's zitterbewegung, for 
graphene stripes of size L in the presence of a transverse mag- 
netic field, B. We show that the interplay between size effects 
and magnetic field can be studied in the continuum limit using 
the Dirac equation coupled to a vector potential. We present 
arguments that show that the spectrum of the problem shows a 
crossover from magnetic field dominated to confinement dom- 
inated as a function of magnetic field or system size. The 
crossover occurs when the radius of the Landau level becomes 
of the order of the width of the system. In the crossover the 



spectrum changes from a/B to linear in B and that the Landau 
plots, that can be measured in a SdH experiment, change from 
dramatically in the presence of a finite system. Our results are 
in excellent agreement with experiments. 
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